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Abstract 

An experimental study of bifurcations associated with stabihty of stationary 

points (SP's) in a parametricahy forced magnetic pendulum and a comparison 

of its results with numerical results are presented. The critical values for 

which the SP's lose or gain their stability are experimentally measured by 

varying the two parameters 17 (the normalized natural frequency) and A (the 

normalized driving amplitude). It is observed that, when the amplitude A 

exceeds a critical value, the normal SP with = (^ is the angle between 

the permanent magnet and the magnetic field) becomes unstable either by a 

period-doubling bifurcation or by a symmetry-breaking pitchfork bifurcation, 

depending on the values of Vt. However, in contrast with the normal SP the 

inverted SP with = vr is observed to become stable as A is increased above 

a critical value by a pitchfork bifurcation, but it also destabilizes for a higher 

critical value of A by a period-doubling bifurcation. All of these experimental 

results agree well with numerical results obtained using the Floquet theory. 
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I. INTRODUCTION 



We consider a permanent magnet of dipole moment m placed in a spatially uniform 
magnetic field B fl]-^. Its motion can be described by a second-order nonautonomous 
ordinary differential equation, 

ie + be + m{BDc + Bac sin ujt) sin ^ = 0, (1) 

where the overdot denotes the differentiation with respect to time, 6 is the angle between 
the magnet and the magnetic field, I is the momemt of inertia about a rotation axis, b is 
the damping parameter, Bdc is the steady DC componet of B, and Bac and u are the 
amplitude and frequency of the sinusoidally time- varying AC component of B, respectively. 
Making the normalization ut —>■ 27rt and 6 —>■ 2ttx, we have 

X + 27r7x + 27r(ri^ + A sin 27rt) sin 2tcx = 0, (2) 

where uq = -y/^f^, ^ = ^> 7 = and A = ^^0^. Note that this is just a normahzed 
equation of motion for the parametrically forced gravitational pendulum (CP) with a verti- 
cally oscillating support 0-0. Hence this magnetic system can be taken as a model of the 
parametrically forced pendulum equation. Hereafter we will call this magnetic oscillator a 
parametrically forced magnetic pendulum (MP). 

The parametrically forced pendulum, albeit simple looking, shows a richness in its dy- 
namical behaviors. As the normalized driving amplitude A is increased, transitions from 
periodic attractors to chaotic attracors and vice versa, coexistence of different attractors, 
transient chaos, multiple period-doubling transitions to chaos, and so on have been nu- 
merically found 1^-0. Some of them have also been observed in real experiments on the 
parametrically forced MP and CP However, so far only the case of = (i.e., 

the case Bj^c = 0) has been studied in the experiments on the MP. 

In this paper, we study experimentally the bifurcations associated with stability of the 
SP's in the parametrically forced MP by varying the normalized natural frequency Q and 
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the normalized driving amplitude A, and then compare the experimental results with the 
numerical results obtained using the Floquet theory |T0|. We first note that the MP has two 



SP's. One is the "normal" SP with {x,x) = (0,0), and the other one is the "inverted" SP 
with (x, x) = (|,0). For the case of the "unforced" simple MP (with A = 0), the normal 
SP is stable, while the inverted SP is unstable. However, as A is increased above a critical 
value, the normal SP loses its stability. Using the Floquet theory, bifurcations occurring 
at such critical values have been numerically studied by one of us (Kim) and Lee 0]. In 
contrast to the normal SP, the inverted SP gains its stability when A exceeds a critical value 



llHl5|. We also study numerically bifurcations associated with stability of the inverted SP 
using the Floquet theory for a comparison with the experimental results. 

This paper is organized as follows. We first explain the experimental setup for the 
parametrically forced MP in Sec. |I|. Bifurcations of the two SP's are then experimentally 
investigated in Sec. |T|. It is observed that as A is increased beyond a critical value, the 
normal SP loses its stability either by a period-doubling bifurcation or by a symmetry- 
breaking pitchfork bifurcation, depending on the values of Q. For the case of the period- 
doubling bifurcation, a new stable symmetric orbit with period 2 is born, while for the case 
of the symmetry-breaking pitchfork bifurcation a conjugate pair of new stable asymmetric 
orbits with period 1 appears. In contrast with the normal SP, the inverted SP is observed to 
gain its stability when A exceeds a 1st critical value A^ by a pitchfork bifurcation, but it also 
destabilizes for a higher 2nd critical value A2 of A by a period-doubling bifurcation. Thus 
the inverted SP becomes stable in the interval between A^ and Ag. Our experimental data 



shows a good agreement between the experimental and numerical results. Finally, Sec. |V 
gives a summary. 

II. EXPERIMENTAL SETUP 

An exploded view of the experimental apparatus is shown in Fig. |1]. The physical MP 
consists of a permanent bar magnet, an aluminum damping plate, and a code wheel, which 
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are coaxially attached to a rotation axis. This MP shows rich dynamical behaviors in a 
spatially uniform magnetic field B, generated by Helmholtz coils. Each component of the 
apparatus is explained in some details below. 

A permanent bar magnet, glued to a rotation axis guided by a small ball in the lowest 
part and by a tiny pin in the highest part, is placed in the center of two sets of Helmholtz 
coils, producing the magnetic field B perpendicular to the rotation axis. The number of 
turns and the radius R for a large set of Helmholtz coils are N = 130 and R = 10.8 
cm, while they are = 144 and R = 5.8 cm for a smaller set of Helmholtz coils. The 
large set of Helmholtz coils is given a direct current Idc to supply a steady DC component 
Bdc of B. On the other hand, the smaller set of Helmholtz coils, which is nested inside the 
large set and driven by a Pasco Model CI-6552A AC power amplifier, provides a sinusoidally 
time- varying AC component of B with amplitude B^c and frequecny u. This amplifier can 
handle currents up to 1 A. Since higher currents are necessary in some range of Q, we also 
use a current booster to increase the current up to 2 A. 

We take into account the effect of the normal component -B^.n of the Earth's magnetic 
field Be perpendicular to the rotation axis, and align both the DC and AC components of 
the applied magnetic field B parallel to Be.h- The angle of the permanent bar magnet is 
experimentally measured from this aligned direction (i.e., the direction of Be-o)- Since the 
effective DC component B^jf of the total magnetic field is given by _Be// = B^c + -B^.n, the 



natural frequency uq in Eq. (|]) becomes ujq = ^yi^e//, where m/I = 2.542 gauss ^ s ^ and 
BE,n = 0.223 gauss 1T|. 



A damping force proportional to the angular velocity can be supplied by an eddy current 
brake. Such eddy-current damping is adjusted by controlling the separation between the 
aluminum damping plate and the horseshoe magnet with a micrometer screw. It is also 
possible to determine a damping parameter k (= 6/2/) and the natural frequency ujq by 
fitting a sampled time-series 6{t) for the unforced case of Bac = to an equation, 6{t) = 



Ce cos {y uJq — k'^ t + 6) {C and 5 are some constants), which is just the angle in the case 
of the underdamped motion of the unforced MP for small angular displacements. 



For data acquisition and experimental control, we use commercial products, "Rotary 
Motion Sensor" and "Signal Interface," manufactured by the Pasco Scientific. A Pasco 
Model CI-6538 rotary motion sensor set consists of a code wheel with 1440 slots (i.e., a 
disk with angular code in the form of sectors which are pervious or impervious to light) 
and an encoder module containing a light-emitting diode (LED) and two photodiodes with 
signal-processing circuitry. As the code wheel moves, the light signal emitted from the LED 
is interrupted by the slots and electrically encoded. Thus a code wheel with 1440 slots can 
generate raw data with a resolution of AO = 0.004 rad. A personal computer equipped with 
a Pasco Model CI-6560 signal interface unit analyses the signals from the rotary motion 
sensor and provides an easy-to-use data set of {6, 6) at a chosen sampling rate. It also 
performs a convenient experimental control. 

III. BIFURCATIONS OF THE NORMAL AND INVERTED SP'S 

In this section we first analyse the bifurcations associated with stability of the SP's in 
the parametrically forced MP, using the Floquet theory. The experimental results for the 
cases of the normal and inverted SP's are then presented and compared with the numerical 
results. 



The normalized second-order ordinary differential equation (|^) is reduced to two first- 
order ordinary differential equations: 



A. Bifurcation analysis based on the Floquet theory 



x = y, 



(3a) 



y = —2'K'yy — 27r(fiQ + Asm2'Kt) sin2'7rx. 



(3b) 



These equations have an inversion symmetry S, since the transformation 



5* : X 



-X, y 



(4) 
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leaves Eq. (|^) invariant. If an orbit z{t) [= {x{t),y{t))] is invariant under S, it is called a 
symmetric orbit. Otherwise, it is called an asymmetric orbit and has its "conjugate" orbit 
Sz{t) [= {-x{t),-ym. 

The Poincare maps of an initial point zq[= {x{0),y{0))] can be computed by sampling the 
orbit points Zm at the discrete time t = m, where m = 1, 2, 3, . . . . We call the transformation 
Zm Zm+i the Poincare (time-1) map, and write z^+i = P{zm)- Note that the normal and 
inverted SP's of the parametrically forced MP, denoted by z^ [= (0,0)] and zi [= (|,0)], 
respectively, are the fixed points (FP's) of the Poincare map P [i.e., P{z) = z for {z = zn, zj)] 
with period 1. 

Here we investigate bifurcations associated with stability of the two normal and inverted 
FP's of P. The linear stability of an FP is determined from the linearized-map matrix DP 
of P at z. Using the Floquet theory [0, we obtain the matrix DP by integrating the 
linearized differential equations for small perturbations as follows. Consider an infinitesimal 
perturbation [6x{t), 6y{t)] to an SP. Linearizing Eq. about the SP, we obtain 

(5) 



(Sx] 




Is \ 

ox 




= j{t) 











where 



J(t) 



(6) 



1 

-Atx'^{VL^ + Asm2-nt)cois,2Txx -27r7 

Here x = and | for the normal and inverted SP's, respectively. 

Note that J is a 2 x 2 1-periodic matrix [i.e., J{t + 1) = J{t)]. Let W{t) = [w^{t), w'^{t)] 
be a fundamental solution matrix with W{0) = I. Here w^{t) and w'^{t) are two independent 
solutions expressed in column vector forms, and I is the 2x2 unit matrix. Then a general 
solution of the 1-periodic system has the following form 



^ 6x{t) ^ 



Sy{t) 



W{t) 



I 6x{0) ^ 



6y{0) 



(7) 



Substitution of Eq. (|^) into Eq. (|^) leads to an initial- value problem in determining W{t): 

6 



W{t) = Jit)W{t), W{0) = I. (8) 

It is clear from Eq. (|^ that VI^(l) is just the hnearized-map matrix DP. Hence the matrix 
DP can be calculated through integration of Eq. (|D over the period 1. 
The characteristic equation of the linearized-map matrix M(= DP) is 

- trM A + det M = 0, (9) 

where trM and detM denote the trace and determinant of M, respectively. The eigenvalues, 
Ai and A2, of M are called the Floquet (stability) multipliers of the PP. As shown in ||17|| , 
det M is calculated from a formula 

detM = e^o''J^\ (10) 

Substituting the trace of J (i.e., trJ = —27ij) into Eq. (pUj), we obtain 

det M = e-^^\ (11) 

Note that M is just a 2 x 2 matrix with a constant Jacobian determinant (less than unity). 
Hence the pair of Ploquet multipliers of an PP lies either on the circle of radius e~'"'^ or on 
the real axis in the complex plane. The PP is stable only when both Ploquet multipliers lie 
inside the unit circle. We first note that they never cross the unit circle except at the real 
axis (i.e., they never have complex values with moduli larger than unity), and hence Hopf 
bifurcations do not occur. Consequently, the PP can lose its stability only when a Ploquet 
multiplier A decreases (increases) through —1 (1) on the real axis. 

When a Ploquet multiplier A decreases through —1, the PP loses its stability via period- 
doubling bifurcation, which leads to the birth of a new stable symmetric orbit with period 
2. On the other hand, when a Ploquet multiplier A increases through 1, it becomes unstable 
via pitchfork bifurcation, which results in the birth of a conjugate pair of new stable asym- 
metric orbits with period 1. Since the newly-born orbits are asymmetric ones, the pitchfork 
bifurcation is also called a symmetry-breaking bifurcation. Por more details on bifurcations. 



refer to Ref. 18 



The stability boundaries of the normal and inverted SP's in some ranges of the Q-A plane 
are determined through numerical calculations of their Floquet multipliers A's. The absolute 
value of A at such stability boundaries is one (i.e., |A| = 1). If A = —1, then the boundary is a 
period-doubling bifurcation line. Otherwise, it is a symmetry-breaking pitchfork bifurcation 
line. We also obtain numerically the bifurcation diagrams and the phase-flow and Poincare- 
map plots at some chosen parameter values for clear visual representation of the bifurcations. 
All of these numerical results are given in the next two subsections (see Figs. for a 
comparison with the experimental results. 

B. Experimental results for the case of the normal SP 

In all the experiments for the normal and inverted SP's, we fix the the driving frequency 
uj in Eq. (^ and the normalized damping parameter 7 in Eq. (Q) as = 2tt and 7 = 0.1, 
respectively. We then control the normalized natural frequency and driving amplitude, fl 
and A, in Eq. (Q) by varying Bdc and Bac, respectively and study the bifurcations of the 
two SP's. 

We consider two ranges of Q for the normal SP, / : 0.2 < Q < 0.5 and // : 0.8 < < 
1.025. For each chosen value of Q, we increase the amplitude A and observe whether the 
SP is stable or not. In order to experimentally determine the stability of the SP, we release 
the MP from rest at a small initial angle displaced from the SP. If the SP is stable, then 
the subsequent motion damps toward the SP. Otherwise, it deviates from the SP. Thus we 
measure experimentally a critical value ^4*^,^ of A, above which the SP is unstable. 

We first study the bifurcations of the normal SP in the 1st range / of Q. As an example, 
consider the case of f2 = 0.4. With increasing A, we carry out the experiments, and measure 
the critical value A^^^. It is observed that for A > A^^^, the SP becomes unstable through 
a period-doubling bifurcation, giving rise to the birth of a new stable symmetric period- 
doubled orbit. For visual representation of the bifurcation, we obtain the bifurcation diagram 
and the phase-flow and Poincare-map plots below. 
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For phase representation, we acquire a data set of [9{t), 9{t)] at a fixed sampling rate 20 
Hz and convert it into a normalized set of {x,y). This whole set of the data is used for a 
phase-flow plot, while a partial set of the data chosen at integral multiples of the external 
driving period T (= 27r) [i.e., t = nT {n = 0, 1,2, ...)] is used for a Poincare-map plot and 
for a bifurcation diagram. 

The bifurcation diagram for Q = 0.4 is shown in Fig. 0(a). The data obtained through 
numerical calculations are also given for a comaprison with the experimental results. Note 
that for the bifurcation diagram, the experimental data represented by the solid circles agree 
well with the numerical data denoted by the solid lines. For reference, the critical values 
obtained through experiments and numerical calculations are ^4*^^ = 0.215 and A*^ = 
0.198 131 ■ ■ ■, respectively. As A is increased above the critical value A*, the normal SP loses 
its stability via period-doubling bifurcation, and a new stable symmetric period-doubled 
orbit appears. Figure ^(b) shows the phase-flow and Poincare-map plots of the symmetric 
period-doubled orbit for A = 0.23. The experimental data for the phase flow are represented 
by the small solid circles, while the two larger solid circles denote the experimental data 
for the Poincare map. These experimental data are also in a good agreement with the 
numerically-computed data for the phase flow represented by the solid line and for the 
Poincare map denoted by the two large circles. 

We also perform the above experiments for many other values of Q in the 1st range 
I and thus measure the critical values A*^p's. The stability diagram of the normal SP 
is shown in Fig. ^. The experimental data for A*^^ are represented by the solid circles, 
and they seem to lie on a smooth stability boundary curve. For a comparison with the 
experimental data, the stability boundary of the SP numerically calculated using the Floquet 
theory is also denoted by the solid line in Fig. ^ This stability boundary is just a period- 
doubling bifurcation line at which a Floquet multiplier of the SP is A = — 1. The period- 
doubling bifurcation line determined through numerical computations lies a little below that 
experimentally determined. That is, the value of A* is somewhat higher than that of 
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This is what one would expect, because in real experiments there exists a frictional force due 
to a contact between the rotation axis and its guiders (ball and pin). As previously noted 
PJ^, one of the main effects of this frictional force is to make the origin of the phase plane 
(i.e., the normal SP) stable up to higher values of the external driving amplitude than the 
numerically-calculated critical value A*^. 

We now study the bifurcations of the normal SP in the 2nd range // of Q (i.e., 0.8 < ^ < 
1.025). As in the above 1st range I of Q, we increase the amplitude A and measure a critical 
value Al^p, beyond which the SP becomes unstable, by releasing the MP from rest at a small 
initial angle displaced from the SP. However, in contrast to the 1st range of Q the normal 
SP is observed to lose its stability through a symmetry-breaking pitchfork bifurcation for 
A = Al^p, which leads to the birth of a conjugate pair of new stable asymmetric orbits with 
period 1. 

As an example, consider the case oi Q = 0.95. The bifurcation diagram for this case 
is shown in Fig. ^(a). The normal SP denoted by the solid circles is observed to become 
unstable through a symmetry-breaking pitchfork bifurcation for A^^^ = 1.25. For A > A^^^, 
a pair of stable asymmetric orbits of period 1 appears. One is represented by the solid circles, 
while its conjugate orbit is denoted by the open circles. Figure ^(b) shows the phase-flow and 
Poincare-map plots of a conjugate pair of symmetry-broken orbits of period 1 for A = 1.4. 
The small solid circles denote the phase flow of an asymmetric "heart-shaped" orbit, while 
the small open circles represent the phase flow of its conjugate "inverted heart-shaped" orbit. 
The data for the Poincare maps of the two symmetry-broken orbits are also denoted by the 
larger solid and open circles, respectively. This symmery-broken case is in contrast to the 
symmetry-preserved case [see Fig. 0(b)] in the 1st range of Q. For a comparison with the 
experimental results, the data obtained by numerical computations are also given in Fig. ^. 
As A is increased above a critical value A*^ (= 1.174 209- ■ ■), the normal SP denoted by 
the solid line becomes unstable through a symmetry-breaking pitchfork bifurcation, giving 
rise to the birth of a conjugate pair of symmetry-broken orbits of period 1. One asymmetric 
orbit is represented by the solid line, while the other one is denoted by the dotted line. 
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As in the experimental case, the data for the Poincare maps of the two asymmetric orbits 
are denoted by the large solid and open circles, respectively. All of these experimental and 
numerical results seem to agree well. 

We also measure the critical values A^^^s for many other values of Q in the 2nd range 
//. Figure |^ shows the stability diagram for this case. The experimental data for A^^^ are 
denoted by the solid circles, while the stability boundary numerically computed using the 
Floquet theory is represented by the solid line. In contrast to the 1st range of Q, the stability 
boundary is a symmetry-breaking pitchfork bifurcation line at which a Floquet multiplier 
of the SP is A = 1. We also note that as in the case of the 1st range /, the value of Agxp 
somewhat higher than that of A*^ because of the frictional force between the rotation axis 
and its guiders. 

C. Experimental results for the case of the inverted SP 

In this subsection, we study the bifurcations associated with stability of the inverted SP 
by increasing A in a range of 0.2 < < 0.5. In contrast to the normal SP, the inverted SP 
is observed to gain its stability when a 1st critical value A* of A is exceeded by a subcritical 
pitchfork bifurcation. However, as A is further increased, the stabilized inverted SP is also 
observed to lose its stability for a 2nd critical value A2 of A through a period-doubling 
bifurcation. Thus the inverted SP becomes stable in the interval between A^ and Ag. 

As an example, we consider the case of i7 = 0.2. The bifurcation diagram for this case 
is shown in Fig. ^(a). The unstable inverted SP denoted by the open circles is observed 
to become stable when A is increased above a 1st critical value A^^^ ^ (= 0.39). Using the 
Floquet theory, the unstable inverted SP denoted by the dashed line is also numerically found 
to gain its stability for A > y4*„ ^ (= 0.289 108 ■ ■ ■) by a subcritical pitchfork bifurcation, 
giving rise to the birth of a conjugate pair of unstable asymmetric orbits with period 1, 
denoted by the dashed lines. However, unfortunately the two symmetry-broken orbits born 
for this subcritical case cannot be experimentally observed, because they are unstable ones. 
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This is in contrast to the supercritical bifurcations occurring for the normal SP in the 2nd 
range II of VL (for a supercritical case, a pair of stable asymmetric orbits is born, which 
can be experimentally observed as shown in Fig. ^). As A is further increased from A*^^^ j^, 
the stabilized inverted SP denoted by the solid circles is observed to lose its stability by a 
period-doubling bifurcation when a second critical value ^exp,2 0.567) is exceeded. For 
A > Al^p2j ^ stable symmetric "butterfly-shaped" orbit of period 2 appears. Small solid 
circles and the two larger solid circles in Fig. P(b) represent the phase flow and Poincare 
map of the symmetric orbit of period 2 for A = 0.61, respectively. It is also numerically 
found that the stabilized inverted SP denoted by the solid line becomes unstable for a second 
critical value A*„2 (— 0.529 159 •■■) through a period-doubling bifurcation, giving rise to 
the birth of a symmetric orbit of period 2 denoted by the solid line. 

We also carry out the above experiments for many other values of Q, and thus measure 
the 1st and 2nd critical values, y4*^p^^'s and Al^p 2^- Figure |^ shows the stability diagram 
of the inverted SP. The experimental data for A^^^^-^ and Al^p 2 represented by the open 
and solid circles, respectively. The inverted SP is observed to become stable in the interval 
between A*^^,^^ and Al^p2- Note also that the width of this interval becomes smaller as 
Q is increased. Hence the stabilization of the inverted SP can be more easily observed 
for small values of Q, compared to the cases of high f2-values. For a comparison with 
the experimental results, numerical data obtained using the Floquet theory are also given 
in Fig. 1^. The lower stability boundary y4*„ ^ denoted by the dashed line is a subcritical 
pitchfork bifurcation line, while the upper stability boundary A^^ 2 denoted by the solid line 
is a period-doubling bifurcation line. We note that the agreement between the experimental 
and numerical results becomes better as Q is decreased. 

IV. SUMMARY 

Bifurcations of the normal and inverted SP's in the parametrically forced MP are experi- 
mentally studied by varying the two parameters Q and A. As A is increased above a critical 
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value, the normal SP is observed to become unstable either by a period-doubling bifurcation 
or by a symmetry-breaking pitchfork bifurcation, depending on the values of In the 1st 
range / of f2 (i.e., 0.2 < fl < 0.5), a new stable symmetric orbit with period 2 appears via 
period-doubling bifurcation, while a conjugate pair of new stable asymmetric orbits with 
period 1 is born via symmetry- breaking pitchfork bifurcation in the 2nd range // of Q (i.e., 
0.8 < fl < 1.025). However, in contrast to this normal SP, the inverted SP is observed to 
become stable when A is increased above a 1st critical value Al by a subcritical pitchfork 
bifurcation. Unfortunately a pair of asymmetric orbits of period 1 born for this subcritical 
case cannot be experimentally observed, because they are unstable ones. As A is further 
increased, the stabilized inverted SP is also observed to lose its stability for a 2nd critical 
value A2 by a period-doubling bifurcation, giving rise to the birth of a new stable symmetric 
period-doubled orbit. Thus the inverted SP becomes stable in the interval between Al and 
A2. When all of these experimental results for the two SP's are compared with the numerical 
results obtained using the Floquet theory, they seem to agree well. 

ACKNOWLEDGMENTS 

This work was supported by the Basic Science Research Institute Program, Ministry of 
Education, Korea, Project No. BSRI-96-2401. 



13 



REFERENCES 

[1] V. Croquette and C. Poitou, J. Phys. Lett. 42, 537 (1981). 
[2] H. Meissner and G. Schmidt, Am. J. Phys. 54, 800 (1986). 
[3] K. Briggs, Am. J. Phys. 55, 1083 (1987). 

[4] V.I. Arnold, Mathematical Methods of Classical Mechanics ( Springer- Verlag, New York, 
1978), p. 113. 

[5] J.B. McLaughhn, J. Stat. Phys. 24, 375 (1981). 

[6] R.W. Leven and B.P. Koch, Phys. Lett. A 86, 71 (1981). 

[7] S.-Y. Kim and K. Lee, Phys. Rev. E 53, 1579 (1996). 

[8] B.P. Koch, R.W. Leven, B. Pompe and C. Wilke, Phys. Lett. A 96, 219 (1983). 

[9] R.W. Leven, B. Pompe, C. Wilke and B.P. Koch, Physica D 16, 371 (1985) 

[10] J. Gukenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifur- 
cations of Vector Fields (Springer- Verlag, New York, 1983), p. 24. 

[11] J.J. Stoker, Nonlinear Vibrations in Mechanical and Electrical Systems (Interscience 
Publishers, Inc., New York, 1950), p. 189. 

[12] H.C. Corben and P. Stehle, Classical Mechanics (Wiley, New York, 1960), 2nd ed., p. 66. 

[13] P.L. Kapitza, in Collected Papers by P.L. Kapitza, edited by D. Ter Haar (Pergamon, 
London, 1965) Vol.2, p. 714. 

[14] V.I. Arnold, Mathematical Methods of Classical Mechanics (Springer- Verlag, New York, 
1978), p. 121. 

[15] J. A. Blackburn, H.J.T. Smith, and N. Gr^nbech- Jensen, Am. J. Phys. 60, 903 (1992); 
H.J.T. Smith and J.A. Blackburn, Am. J. Phys. 60, 909 (1992). 



14 



[16] By decreasing the value of Boc from some positive value, we first get a data set of 
uq^Bdc) for several values of -B_dc through the fitting procedure described in the 4th 
paragraph in Sec. 0. We can then obtain the values of BE,n and m/I by fitting the 
experimental data set of ujo{Bdc) to the equation, ujq = \J^{Bdc + Be^u)- 

[17] S. Lefschetz, Differential Equations: Geometric Theory (Dover Publications, Inc., New 
York, 1977), p. 60. 

[18] J. Gukenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifur- 
cations of Vector Fields ( Springer- Verlag, New York, 1983), Sec. 3.5. 



15 



FIGURES 

FIG. 1. Exploded view of a parametrically forced MP. Labeled components are the aluminum 

damping plate (^4), the horseshoe magnet controlled by a micrometer (B), the code wheel (C), the 
encoder module (D), the Helmholtz coil for production of a steady DC component of a spatially 
uniform magnetic field B (E) , the Helmholtz coil for production of a sinusoidally time- varying AC 
component of B (F) and the permanent bar magnet (G). 

FIG. 2. Period-doubling bifurcation of the normal SP for Jl = 0.4. As shown in the bifurcation 
diagram (a), the normal SP denoted by the solid circles becomes unstable for a critical value A^^^ 
(= 0.215) by a period-doubling bifurcation, which results in the birth of a symmetric period-doubled 
orbit denoted by the solid circles. The experimental data for the phase flow for A = 0.23 are denoted 
by the small solid circles in (b). Numerical data denoted by the solid lines are also given in both 
(a) and (b). The data for the Poincare maps in (b) are represented by the two large solid circles 
for both the experimental and numerical cases. For more details see the text. 

FIG. 3. Stability diagram of the normal SP in the 1st range I oi Q (i.e., 0.2 < O < 0.5). 
The experimental data for the critical values A*^pS, above which the SP becomes unstable, are 
represented by the solid circles. The stability boundary numerically determined using the Floquet 
theory is denoted by the solid line, and it is just a period-doubling bifurcation line. 
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FIG. 4. Symmetry-breaking pitchfork bifurcation of the normal SP for U = 0.95. For a critical 
value Al^p (= 1.25), the normal SP denoted by the solid circles becomes unstable by a symme- 
try-breaking pitchfork bifurcation, as shown in the bifurcation diagram (a). As a result of the 
bifurcation, a conjugate pair of asymmetric orbits of period 1 appears; one is denoted by the solid 
circles, while the other one is represented by the open circles. The experimental data for the phase 
flow of the two symmetry-broken orbits for A = lA are denoted by the small solid and open circles 
in (b) , respectively. Numerical data for the two symmetry-broken orbits are also given in both (a) 
and (b); one is denoted by the solid line, while the other one is represented by the dotted line. The 
data for the Poincare maps in (b) are represented by the large solid and open circles for both the 
experimental and numerical cases. For more details see the text. 

FIG. 5. Stability diagram of the normal SP in the 2nd range // of fl (i.e., 0.8 < Q < 1.025). 
The experimental data for the critical values A^^p's, above which the SP becomes unstable, are 
represented by the solid circles. The stability boundary numerically determined using the Floquet 
theory is denoted by the solid line, and it is just a symmetry-breaking pitchfork bifurcation line. 

FIG. 6. Bifurcations of the inverted SP for O = 0.2. The bifurcation diagram for Q = 0.2 is 
shown in (a). In contrast to the normal SP, the inverted SP denoted by the open circles is observed 
to become stable as A is increased above a 1st critical value A^^p ^ (= 0.39). However, as A is 
further increased, the stabilized inverted SP denoted by the solid circles is also observed to lose its 
stability for a 2nd critical value ^exp,2 (= 0.567) through a period-doubling bifurcation, giving rise 
to the birth of a symmetric period-doubled orbit denoted by the solid circles. The experimental 
data for the phase flow of the period-doubled orbit for A = 0.61 are denoted by the small solid 
circles in (b). Numerical data denoted by the lines are also given in both (a) and (b); a stable orbit 
is denoted by a solid line, while an unstable orbit is represented by a dashed line. The data for the 
Poincare maps in (b) are represented by the two large solid circles for both the experimental and 
numerical cases. For more details see the text. 
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FIG. 7. Stability diagram of the inverted SP in the range of 0.2 <Q < 0.5. The experimental 
data for the 1st and 2nd critical values, A^^^ -^ and Al^p 2i ^^e represented by the open and solid 
circles, respectively. The lower and upper stability boundaries, numerically computed using the 
Floquet theory and denoted by the dashed and solid lines, are the symmetry-breaking pitchfork 
and period-doubling bifurcation lines, respectively. 
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